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Abstract 

From considering the effect of 7-7 interactions on recently observed TeV 
gamma-ray spectra, improved limits are set to the density of extra- 
galactic infrared (IR) photons which are robust and essentially model- 
independent. The resulting limits are up to two orders of magnitude 
more restrictive than direct observations in the 0.025-0.3 eV regime. 
These limits are used to improve constraints on radiative neutrino decay 
in the mass range above 0.05eV and on Very Massive Objects (VMOs) 
as providing the dark matter needed to explain galaxy rotation curves. 
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1 Introduction 



The extragalactic background IR field potentially contains a wealth of informa- 
tion relevant to both cosmology and particle physics (see for example (![]), but 
has so far eluded any conclusive detection. Direct measurement of the charac- 
teristics of this radiation are frustrated by the dominance of local, galactic IR 
sources. TeV gamma-ray astronomy provides a means to study the IR back- 
ground indirectly, free of such complications, by looking for modifications to 
high energy gamma ray spectra due to interactions with this background field 
via the process 77 — ► e + e~ . This idea was first noted by Gould and Schreder 
[0] and has been more recently restated by Stecker and de Jager |H|. 

Previous IR studies making use of this phenomenon have generally used 
spectral data from the active galactic nucleus (AGN) Mrk 421 taken by the 
EGRET m and Whipple experiments at GeV and TeV energies, respec- 
tively. These studies have relied on extrapolating flux measurements between 
these two energy regimes in an effort to determine the inherent shape of the 
source spectrum 0] HI HI- This approach has several difficulties: 1) Since 
the flux level is known to be highly variable at TeV energies, the relative nor- 
malization between EGRET and Whipple flux measurements is not known; 2) 
The flux of IR photons cannot be assumed to be zero in the intervening energy 
range; and 3) The shape of the inherent source spectrum is not known and 
cannot be assumed to necessarily follow a single power law from GeV to TeV 
energies. Indeed, current best fit values for Mrk 421 at TeV energies indicate 
a steeper spectrum than is reflected by the EGRET data [0, suggesting that 
the spectral shape changes in the intervening energy range. In addition, most 
of the analyses have assumed a particular model to describe the shape of the 
IR photon spectrum. Even a more recent attempt at a less mo del- dependent 
limit using data from Mrk 501 nevertheless assumes a power-law source spec- 
trum, an arbitrary IR flux normalization and constant IR energy densities 



over energy ranges spanning up to an order of magnitude in extent [II]. Lim- 
its derived in this manner may not apply if the background IR spectrum is of 
a different shape and, in particular, cannot necessarily be used to constrain 
mechanisms which might produce sharp features in the spectrum on a scale 
smaller than that of the assumed model. This current paper will attempt to 
define more robust limits based on recent high energy observations of AGN 
which do not rely on extrapolation from lower energies or the assumption of a 
strict power-law source spectrum and are independent of the IR model down 
to scales corresponding to a factor of 2 change in IR energy. These limits 
will then be used to constrain significantly the possible contribution of Very 
Massive Objects (VMOs) to the dark matter problem and also to place im- 
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proved bounds on the radiative decay of massive neutrinos in the mass range 
above ~ 0.05 eV, which is in the regime of interest to the atmospheric neutrino 
problem. 



2 Derived Upper Limits to the Background IR 

Two active galaxies, Mrk 421 and Mrk 501, have thus far been seen to produce 
TeV gamma-ray emission by independent experimental groups. These active 
galaxies are both of the blazar class and are located at redshifts of 0.031 and 
0.033, respectively. Both objects have exhibited significant flux variability at 
TeV energies, changing by an order of magnitude or more with characteristic 
timescales of less than an hour (in the case of Mrk 421, as short as 15 minutes). 
Recent analyses of data from Mrk 421 by the Whipple group || (|TIJ indicate a 
differential spectrum with an ~ E~ 2M power law extending to energies beyond 
5 TeV, possibly as high as 8 TeV. Similarly, the HEGRA collaboration has 
measured the "high-state" spectrum of Mrk 501 [0] and has found it to be 
described by an ~ E~ 2 - 5 power law extending to energies as high as 10 TeV. 

For the purposes of this paper, the available data on active galaxies at a 
redshift of ~ 0.03 will thus be taken to indicate a spectral shape which is 
bounded at the 90% confidence level between the power laws E~ 2 ' 3 and E~ 2 ' 8 
and normalized by an arbitrary flux constant. It will be assumed that the 
measured spectra continue up to energies of 10 TeV. However, due to possible 
systematic uncertainties in energy scale of ~30%, more conservative limits will 
also be computed assuming a maximum gamma-ray energy of 6 TeV. Finally, 
it will be assumed that the inherent differential source spectrum is steeper 
than an E~ 2 power law, although the specific shape need not be a power law 
itself. The limits derived are relatively insensitive to this latter choice. For 
example selecting E~ x ^ as a bound was found to change the inferred limits by 
less than 30%. 

For ease of computation, the high energy spectrum is taken to consist of 
differential flux measurements at the following discrete energies: 0.4, 0.75, 
1.0, 1.5, 2.0, 3.0, 4.0, 6.0 and 10.0 TeV. The spectrum of IR photons is ap- 
proximated by a series of separate intervals spanning the energies 0.025-0.042, 
0.042-0.083, 0.083-0.167, 0.167-0.33, 0.33-0.625 and 0.625-1.0 eV. Within each 
interval, the value of e 2 n(e) is taken as constant, where s is the energy of the 
IR photon and n is the number density of photons (i.e. photons/cm 3 /eV). 
The value of e 2 n(e) among these intervals can thus be varied to reflect differ- 
ent background IR models, and the corresponding effect of 7-7 absorption on 
TeV gamma-rays calculated according to the prescription described by previ- 
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ous papers 0,0. For these calculations, a conservative value for the Hubble 
constant of 85 km/s/Mpc is assumed [13 . 

To calculate upper limits then, the inherent AGN differential spectrum 
was allowed to assume any shape steeper than E~ 2 and the values of e 2 n{e) 
among IR intervals were allowed to take on any values subject to the following 
constraints: 1) that these values are lower than the directly determined upper 
bounds imposed by other experiments (shown in figure 1); and 2) that the 
attenuation implied by the IR model would not cause the shape of the observed 
TeV spectrum to deviate outside of the power-law bounds previously defined 
based on the data. This approach essentially anchors the lower energy TeV 
data to the appropriate observational upper bounds and then extends these 
bounds based on the shape and extent of the AGN spectra at higher energies. 
From these considerations, the maximum value of energy density permitted in 
each IR interval was determined. The resulting limits are summarized in table 
1 and are also plotted in figure 1 for assumed maximum gamma-ray energies of 
6 TeV and 10 TeV. Note that this choice of maximum energy primarily effects 
only the two lowest energy IR intervals. 





Upper Limit to e 2 n(e) (eV/cm 3 ) 


Energy Interval (eV) 


E max = 10 TeV 


Emax = 6 TeV 


0.025-0.042 


0.0064 




0.042-0.083 


0.0027 


0.094 


0.083-0.167 


0.0027 


0.0038 


0.167-0.33 


0.0035 


0.0037 


0.33-0.625 


0.0052 


0.0052 


0.625-1.0 


0.012 


0.012 



Table 1: Derived upper limits to e 2 n(e) 



In the range of 0.025 — 0.33 eV, the limits presented here are up to 2 orders 
of magnitude more restrictive than those of DIRBE Note that these limits 



apply only to models which would produce features in the IR spectrum on a 
scale comparable to or larger than the interval widths used in the calculation. 
Monoenergetic line features, for example, are not necessarily constrained by 
these limits. 
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3 Constraints on VMO Abundance 



Under certain cosmological scenarios, density fluctuations in the early universe 
at a redshift of ~1000 could give rise to the pregalactic formation of massive 
stars (200 < M/M < 10 5 ) at redshifts in the range of ~ 100-300 pf. Such 
Very Massive Objects (VMOs) would collapse to black holes after their main- 
sequence phase without any significant metal ejection. It has been suggested 
that the remnants of VMOs, which would ultimately cluster with galaxies, 
might provide an explanation for the dark matter associated with galaxy ro- 
tation curves fl6| , |I7| . Such rotation curves indicate that the amount of dark 
matter is in excess of 5 percent of the critical density, while other studies of 
galaxy motion and hot gas in clusters of galaxies suggest a value as high as 20 
percent jl3 l . 



While invisible now, VMOs would have produced a significant flux of IR 
photons during their burning phase. The peak number density of these photons 
is given pj by: 



e 2 pk n v (e pk ) ~ 40 



ttyh 2 
l + z v 



30 

where e p k is the IR photon energy (in eV) corresponding to the peak of 
the VMO contribution, ny is the photon number density due to VMOs in 
cm _3 eV -1 , fly is the mass-density of VMOs expressed fraction of the 
critical density, zy is the formation redshift of the VMOs, and the Hubble 
constant is given by Ho = 100/i km/s/Mpc. Thus, the IR limits derived pre- 
viously may be used to place upper bounds on fly. These bounds are shown 
in figure 2 as a function of the assumed VMO formation redshift for h—0.65 
(which is conservative for this limit). Note that for either maximum gamma- 
ray energy assumed, the limit on fly is less than 0.05 for the redshift range of 
100-200 and less than 0.065 for the range 200-300. If the maximum gamma-ray 
energy is taken to be 10 TeV, the limit is less than 0.05 for the the entire 100- 
300 range. Hence, assuming minimal reprocessing by dust grains, these limits 
appear to rule out VMO models as providing the explanation for observed 
galaxy rotation curves. 
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4 Constraints on Radiative Neutrino Decay 



Several experiments have observed an apparent descrepency in the ratio of u ( 
to Up fluxes expected from cosmic ray interactions in the atmosphere [18| [19 



[ p0|1 . A possible explanation for this effect may involve 2-component neutrino 



oscillations with Am 2 ~ O.OleV 2 and, thus, a neutrino mass in the regime of 
0.1 eV (provided the mass of the lighter neutrino is comparable to or smaller 
than this value). If neutrinos have mass and different flavor eigenstates may 
mix, then massive neutrinos might decay via the process vjj — > + 7, where 
m UH > m UL . For neutrino masses in the ~eV regime, such a radiative de- 
cay process would contribute to the infrared background. Although mono- 
energetic photons would be initially produced from a 2-body decay, the photon 
spectrum would be broadened for relic neutrinos which decay over a range of 
redshifts such that the background density can be constrained by the limits 
presented previously. 

Some unstable relic neutrinos with a lifetime longer than that of the uni- 
verse would still be decaying now and would have produced a background 
photon density given by: 

nJe) ~ 1.1 x 10 12 -^- 
r v ht 

where e is half the rest mass of the unstable neutrino in eV, n 7 is the photon 
number density due to neutrino decay in cm~ 3 eV~ 1 , B r is the branching ratio 
for the radiative decay mode and r v is the neutrino lifetime in years |21 |. 



Implied lower bounds for r u /B r in the neutrino mass range O.leV^ < m v < leV 
are shown in figure 3 based on the IR limits presented in this paper for h = 0.85. 



Previous constraints in this mass range from the work of Ressell and Turner |21 
are shown for comparison. This represents an improvement of the constraints 
in some cases by more than 2 orders of magnitude. 

Similarly, neutrinos with lifetimes shorter than the age of the universe but 
longer than the age at which decoupling of matter and radiation occurred will 
result in a background photon density given by: 



n 7 (e) ~ 115 — 



where, in this case, 



-m„c 2 \ / T„\/TT h \ 2 ^" 



2 



eV J \ 6.5 x 10V / 



6 



where Q is the mass-density of the universe in units of the critical density. IR 
limits given here have thus been used to define exclusion contours which can 
be conveniently characterised by 



where 0.05eV < m < 2.0eV for an assumed maximum gamma-ray energy of 
10 TeV, and Q.leV < m < 2.0eV for a maximum energy of 6 TeV. Once again, 
this represents an improvement in branching ratio sensitivity of approximately 
two orders of magnitude. 

5 Discussion 

Recent data from very high energy gamma-ray observations of active galaxies 
has been used to set improved upper limits to the infared background radiation 
which are robust and relatively model-independent. These limits cast doubt 
on VMO models as providing the explanation for dark matter in galactic halos 
and tighten constraints on the radiative decay of relic neutrinos in the mass 
range near 0.1 eV by approximately two orders of magnitude. 

Substantial improvement in the ability to study the infrared background 
via 7-7 interactions is expected with the next generation of instruments which 
will allow the behavior of more distant sources to be examined. In fact, the 
procedure used to derive the limits given in Table 1 may also be used to 
determine lower limits to the maximum distance out to which TeV gamma-ray 
telescopes can be expected to probe. This has been done for the gamma-ray 
energy range of 0.4-10 TeV, assuming that a source remains "visible" out to 
an optical depth of 2 and that the majority of the IR was produced prior 
to the epoch of the sources under study. IR and UV upper bounds from 
direct observations in the energy range greater than 2 eV have also been used 
to calculate this minimum distance for gamma-ray energies of 0.05 and 0.1 
TeV. The results are shown in figure 4. The current generation of ground- 
based, gamma-ray telescopes have a typical lower energy threshold of ~0.5 
TeV and are thus expected to be able to see sources possessing redshifts up 
to or beyond z—0.1. The next generation of instruments is expected to have 




for the region defined by 
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an energy threshold in the region of 0.05 to 0.1 TeV, and will therefore be 
able to see out to a redshift of at least z—0.5. This represents a substantially 
larger volume of space than has so far been explored at TeV energies. This, 
and the known preponderance of AGN at larger redshifts, indicates a bright 
future for ground-based, gamma-ray astronomy and strongly suggests a drive 
towards higher duty cycles, larger apertures and greater dynamic ranges for 
such instruments to explore this regime thoroughly. 

This work has been supported in part by PPARC, Forbairt, the US De- 
partment of Energy, NASA and the Smithsonian Institution. 
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Figure 1: Upper limits to the background infrared background: "+" symbols 
are from Matsumoto et al. [p2| ; the diamond is from Dube et al. |23|| : the 
triangle is from Toller | 24}[ ; squares are from Hauser a/. [14]; circles are 



from Kashlinsky et al. [P3J; crosses are from Maucherat-Joubert et al. 



the solid lines show the limits derived in this work for a maximum gamma-ray 
energy of 10 TeV, while the dashed lines correspond to a maximum energy of 
6 TeV. 
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Figure 2: Upper limits to the abundance of Very Massive Objects based on 
bounds to the IR background derived in this work. The limit is expressed as 
a fraction of the critical density and is plotted as a function of the assumed 
formation redshift of the VMOs. The solid line corresponds to an assumed 
maximum gamma-ray energy of 10 TeV and the dashed line to a maximum 
energy of 6 TeV. 
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Figure 3: Lower bounds to the ratio of the neutrino lifetime to the branching 
ratio for radiative decay based on this current work for the case where the 
neutrino lifetime is greater than the current age of the universe. The solid line 
shows the limits derived here assuming a maximum gamma-ray energy of 10 
TeV, while the dashed line coresponds to a maximum energy of 6 TeV. Open 
squares show the limits derived previously by Ressell and Turner. 
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Figure 4: Lower limits to the maximum redshift out to which TeV gamma-ray 
sources are visible as a function of primary gamma-ray energy, based on limits 
to the background density of IR/UV photons. 
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